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^^ ' Abstract. We analyze self-consistent A^-body simulations of the Milky Way disk and the 

r^ ■ ongoing disruption of the Sagittarius dwarf satellite to study the effect of Sagittarius tidal 

Qh! debris on dark matter detection experiments. In agreement with significant previous work, we 

Q ' reiterate that the standard halo model is insufficient to describe the non-Maxwellian velocity 

distribution of the Milky Way halo in our equilibrium halo-only and halo/galaxy models, and 

^ , offer suggestions for correcting for this discrepancy. More importantly, we emphasize that the 

dark matter component of the leading tidal arm of the Sagittarius dwarf is significantly more 

extended than the stellar component of the arm, since the dark matter and stellar streams 

^ . are not necessarily coaxial and may be offset by several kpc at the point at which they impact 

t^^ I the Galactic disk. This suggests that the dark matter component of the Sagittarius debris is 

\^ • likely to have a non-negligible influence on dark matter detection experiments even when the 

^^ . stellar debris is centered several kpc from the solar neighborhood. Relative to models without 

ff^ I an infalling Sagittarius dwarf, the Sagittarius dark matter debris in our models induces an 

^^ ' energy-dependent enhancement of direct search event rates of as much as ~ 20 — 45%, an 

energy-dependent reduction in the amplitude of the annual modulation of the event rate by 

as much as a factor of two, a shift in the phase of the annual modulation by as much as ~ 20 

days, and a shift in the recoil energy at which the modulation reverses phase. These influences 

/\ • of Sagittarius are of general interest in the interpretation of dark matter searches, but may 

C^ , be particularly important in the case of relatively light {rriy, < 20 GeV) dark matter because 

the Sagittarius stream impacts the solar system at high speed compared to the primary halo 

dark matter. 
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1 Introduction 

Numerous experiments aim to detect particle dark matter directly by measuring the rate of 
extremely rare nuclear recoil events from the elastic scattering of weakly-interacting massive 
particles (WIMPs; for context on this process, see [1-3]). To interpret properly any putative 
direct detection signal or the limits on particle properties implied by any null searches, a de- 
tailed understanding of the phase space distribution of dark matter in the solar neighborhood 
is required ([4-6]; see also [7] for general review of the topic). Predictions for the rate of scat- 
tering events in direct search experiments are usually made under standard assumptions of a 
canonical local WIMP density, po ~ 0.3 — 0.4 GeV/cm^, and a WIMP velocity distribution 
described by a Maxwellian function with a three-dimensional dispersion, cjsd ^ 270±70 km/s 
[8]. 

Current experiments probe the parameter space of WIMP mass and scattering cross- 
section addressing many interesting dark matter candidates, including the lightest super- 
partner particles in supersymmetric theories [9]. Tentative signals suggest that the WIMPs 
comprising the cosmological dark matter may be particles with mass on the order of rriy. ~ 
5-20 GeV/c2 (e.g., the CRESST, DAMA/LIBRA, and CoGeNT experiments in [10-12], re- 
spectively). In this mass range, nuclear recoils with sufficient energy to be detectable require 
larger relative velocities than dark matter particle candidates in the higher mass range of 
m^^ ~ 100 — 1000 GeV/ that has been more widely explored in recent years (as demonstrated 
by Figure 1 in ref. [13]). The requirement of larger relative velocities renders scattering rates 
significantly more sensitive to the high-velocity tail of the dark matter velocity distribution 
than has been considered by many recent studies [14, 15]. 

The slight velocity anisotropy and significant spatiotemporal variation of the velocity 
distribution [16-18] in a Milky Way-sized dark matter halo has been shown to affect signifi- 
cantly the predicted direct detection event rate compared to that obtained from a standard 
Maxwellian velocity distribution [19, 20]. The possibility of dark matter streams in the solar 
neighborhood has induced theoretical investigations of the effect of a coherent population of 
WIMPs on the annual modulation of the event rate observed by direct detection experiments 
[21-26]. Meanwhile, high-resolution cosmological numerical simulations of halo formation 
have revealed that there are many kinematically-cold streams comprising the halo as well as 
myriad debris inflows from stripped subhalos [27, 28]. These streams are associated with past 



mergers that built the Galactic mass to date, and are almost all of insufficient mass density 
to affect event rates measurably at the Earth's location in the Milky Way halo [29, 30]. 

Previous numerical studies that have addressed the issue of direct detection and ve- 
locity distributions in high-resolution A^-body simulations (e.g., refs. [18, 19]) have largely 
been limited to analyses of simulations including only dark matter. These simulations do 
not account for the Milky Way disk and its cosmological growth and evolution to the present 
day. Such experiments represent only individual samples of the broad statistical ensemble of 
merger histories that lead to the formation of a Milky Way-sized dark matter halo. Conse- 
quently, these simulated halos lack specific structures known to exist in the Milky Way, such 
as the ongoing Sagittarius accretion and tidal stream, or other features that may have yet 
to be discovered. Therefore, the problem of mapping the results of cosmological numerical 
simulations onto specific predictions for the actual Milky Way halo, especially at the solar 
neighborhood, is a distinct challenge. 

Likewise, cosmological simulations designed to model the growth of a disk galaxy are not 
yet capable of accounting for the fundamental structures in a Milky Way-like galaxy. Such 
simulations generally do not produce thin and dynamically-cold stellar disks analogous to 
that of the Milky Way (as demonstrated by ref. [31] for a variety of computational techniques 
and energy feedback algorithms). Moreover, the computational expense of these campaigns 
imposes limits on the degree to which they can account for constraints on the unique mass 
assembly history of the Milky Way itself [32, 33]. As an example, several hydrodynamical 
simulations have produced a Galaxy-sized host dark matter halo in a cosmological box, in 
concert with a significant, coherently-rotating "dark disk" component due to large accretions 
at late times [34-36] - although these components are widely found in a ACDM universe, 
the thin and cold Milky Way has had a quiescent accretion history that is incompatible with 
such recent minor mergers [37]. 

Modeling the dark matter and luminous components of the Milky Way in equilibrium 
with each other, in an isolated simulation calibrated to specific characteristics of the Milky 
Way (as in the formalism of [38]), complements cosmological simulations in a number of ways. 
Isolated simulations are presently rather well-constrained by observational data on Galactic 
structure and substructure and therefore may represent many aspects of the Milky Way 
more faithfully than an individual sample from a highly-stochastic ensemble of merger and 
accretion histories of Milky Way-sized halos in dark matter simulations that neglect baryons. 
Such isolated calculations represent one method for estimating the infiuences of significant 
baryonic structures, such as the Galactic disk, and controlling for the unique merger history 
of the Milky Way. Isolated simulations also provide a framework for incorporating important 
accidental features such as the Sagittarius dwarf galaxy (Sgr) and its associated tidal stream. 
The disadvantage of this approach is that it does not preserve the cosmological construction 
of the host halo, in which virialized streams associated with high-redshift merger activity 
would be remnants of galaxies with so few stars that their present-day spatial distribution 
cannot be determined; these mergers are not self-consistently treated in an approach explicitly 
designed to model the current structure of the Milky Way, in equilibrium with a smooth dark 
matter host halo. In the absence of good constraints on the late-time accretion history of 
the Galaxy, as determined observationally by mapping surveys, treating the Milky Way halo 
as an isolated and equilibrated disk/bulge/halo model at redshift 2; = is a complementary 
method to high-resolution cosmological simulations. 

In this paper, we analyze a set of isolated simulations designed to represent the Milky 
Way disk along with the Sagittarius galaxy impact on the disk and the associated Sagittarius 
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Figure 1. Upper panels: distributions of stars and dark matter from the disrupting Sagittarius dwarf 
satellite galaxy {green and gray particles, respectively), in our light Sgr and heavy Sgr models. In this 
perspective, the Milky Way disk plane is denoted by concentric blue rings at 5-kpc radial intervals, 
to a Galactoccntric distance of 25 kpc. Lower panels: surface density maps of Sgr stars and dark 
matter through the disk mid-plane for each model, computed in a slice with depth of 2 kpc, similar 
to the width of the Galactic disk. Note that the leading stellar stream does not pass directly through 
the solar neighborhood, although a significant amount of dark matter belonging to that stream is 
found near the Sun (as shown by red particles in the upper inset panels, and by grayscale shading in 
the lower panels). Indeed, the projected contours of stellar density and dark matter density are not 
concentric, particularly in the more heavily-stripped light Sgr model. 



tidal debris. We demonstrate that neglecting important modifications to A^-body predictions, 
such as the presence of a large stellar Galactic disk and the ongoing Sagittarius dwarf merger, 
both structures that are known to exist, may well result in systematic misestimations of the 
expected event rates and annual modulation signatures in direct detection experiments. The 
stellar disk of the Milky Way itself, being in approximate equilibrium with its dark matter 
host halo, has concomitantly drawn the near-midplane region of the dark matter halo into 
a phase-space distribution peaked at higher velocity and with larger deviations from the 
Maxwellian form than found in dark matter-only simulations. These results follow from 
models of disk-halo equilibria formulated in ref. [38] and are related to numerous previous 
statements of uncertainty in the local dark matter density and velocity distribution in Milky 
Way models, as in the reconstructions of [39, 40] among many other efforts. 



Our most interesting results address the influence of the Sagittarius debris on dark 
matter detection experiments. Some observational measurements of the tidally-stripped stars 
suggest that the leading stellar arm of the Sgr tidal stream may fall several kiloparsecs away 
from the solar position in the Galactic plane (e.g., [41, 42]); although recent detections 
of coherently-moving stellar populations in the solar neighborhood seem to rule out large- 
scale flows with Sagittarius-like vertical velocities [43, 44], percent-level streams outside one 
or two kiloparsecs from the Sun are poorly constrained and difficult to disentangle using 
tracer-sampling techniques. Observationally-viable models for the Sgr infall exhibit several 
interesting features relevant to contemporary and future direct search experiments, because 
the models we present here treat the Sgr dwarf in a cosmological context, in the sense that we 
assume that the Sagittarius galaxy is itself embedded in a dark matter halo as it merges with 
the Milky Way. Utilizing simulations of the Sagittarius infall, we emphasize that although 
the primary stellar component of the Sagittarius debris stream may be several kiloparsecs 
away, the associated dark matter stream is always significantly more spatially extended than 
the stellar stream (as shown in Figure 1) and more coherent in velocity space. In fact, 
our models illustrate that the stellar and dark matter streams need not even be spatially 
concentric in realistic models. In other words, the peak of the stellar density of the stream 
may be displaced from the peak in dark matter density associated with the stream by several 
kiloparsecs. Both of these properties of simulated Sagittarius infall models suggest that it 
is not necessary for the solar neighborhood to be closer than several kiloparsecs away from 
the primary stellar stream of Sagittarius debris in order for the stream to affect significantly 
dark matter detection rates. We use the dark matter streams realized in our simulations to 
model the dark matter velocity distribution in the solar neighborhood. Indeed, we find that 
the near coincidence of the solar neighborhood and the debris of the Sagittarius satellite's 
dark matter halo creates a peak in the high-end of the dark matter velocity distribution. 
These Sgr particles can lead to significant effects on the energy dependence of dark matter 
detection rates, as well as the amplitude and phase of the annual modulation of event rates, 
particularly for relatively low-mass dark matter candidates {m^ < 20 GeV). 

Our detailed calculations involve velocity distributions drawn from high-resolution nu- 
merical experiments, involving an isolated Galaxy-analog stellar disk/bulge/halo model and 
self-consistent treatment of the Sagittarius infall. In §2, we discuss the relevant features of 
simulations introduced by [45] and describe the analysis implying that significant alterations 
may be necessary to the predictions of event rates in direct detection experiments. Our re- 
sults are presented in §3, reserving §4 for discussion of their application and potential future 
prospects regarding detectors with low recoil-energy thresholds and more sensitive annual 
modulation constraints. 

2 Methods 

2.1 Simulations of the Galactic Environment 

All analyses in this work are based on simulations described in [45], and we refer the reader 
there for technical details. In brief, we describe four high-resolution experiments based on 
self-consistent halo/galaxy models, using throughout the notation "/lost halo" to refer to a 
standard NFW dark matter halo [46] with a mass profile roughly consistent with that of the 
Milky Way, and ^^halo + disk" for the model with an identical halo generated in equilibrium 
with a Galactic-analogue disk and central bulge with a global density structure consistent 
with that of the observed Milky Way stellar disk [38], and a total-mass surface density 



at the solar location that is constrained to match the observations of ref. [47] and therefore 
consistent with recent refinements to this value [48] . For the purposes of the present paper, in 
the context of phase-space behavior on small scales, we note that this three-component model 
represents an approximate, self-consistent solution to the coupled collisionless Boltzmann and 
Poisson equations, and that all A'-body simulations were performed using many millions of 
particles with a force-softening resolution equal to one parsec. The disk is stable against 
large-scale perturbations, and develops only a weak bar in isolated evolution over a timescale 
of 3 — 5 Gyr. The disk remains globally static in density and velocity structure outside a 
few kiloparsecs from the Galactic center, ensuring that the solar neighborhood in this model 
remains unaffected by artificial baryonic asymmetries related to spiral-arm evolution. 

Here and throughout, we define the solar neighborhood as a wedge in the radial range 
8 < -R0 < 9 kpc and extending ±2.5 (1.5) kpc in the vertical (tangential) direction. Wedges 
of this size are used to define the kinematic properties of the halo and stream dark matter par- 
ticles because a large number of particles is needed to probe velocity structure. The relative 
densities of the dark matter and stream components can be determined within significantly 
smaller volumes. In each case, we have verified that we derive consistent results when we vary 
the size of the wedge by ~ 1 — 2 kpc in the tangential direction and up to ~ 10 kpc in the verti- 
cal direction. The insensitivity of our results follows from the relatively mild spatial variations 
in stream density exhibited in Figure 1. Following the example of ref. [14], the rest-frame on 
Earth is taken to include the galactic rotation (a purely tangential Vc = 220 km/s in the host 
halo model, and calculated according to the stellar disk's rotation in the halo + disk and 
related models), as well as the peculiar solar motion {Uq,Vq,Wq) = (11.1,12.2,7.3) km/s 
according to [49], and finally the Earth's orbital motion around the Sun as prescribed by 
[50]. The position of the impact of the Sagittarius stream on the Galactic disk is somewhat 
uncertain, so there is some arbitrariness in our identification of the solar neighborhood. We 
identify the solar neighborhood such that the Sun has approximately the correct position rel- 
ative to the Sagittarius dwarf remnant, and we return to uncertainty in the relative position 
of the Sun with respect to the stream in § 3.2. 

As in [45], we examine two numerical experiments in which the Galactic disk is impacted 
by cosmologically-realistic progenitors to the Sagittarius dwarf galaxy, resulting in significant 
spirality and ring-like structure in the host Milky Way. The dark matter halo masses of the 
two satellite models roughly bracket the expected range motivated by cosmological abundance 
matching techniques [e.g., 51]: "light Sgr''' assumes a Sagittarius progenitor mass of Afsgr = 
IQ-'^'^'^AfQ, while "heavy Sgr'^ has a progenitor mass of Mggr = IQ^^Mq. This cosmologically- 
plausible range is supported by the dynamical reconstruction of the progenitor as reported 
in ref. [52], in which the mass immediately prior to tidal disruption is estimated to be at 
least ~ IO^^'Mq, a value consistent with the total masses contained within and truncated at 
the two model subhalos' initial Jacobi tidal radii when our simulations begin. As reported 
in [45], both models are consistent with the approximate characteristics of the Sagittarius 
tidal debris as mapped by the Two-Micron All Sky Survey and the Sloan Digital Sky Survey 
among other efforts [53-55]. In particular, these evolved debris distributions reproduce the 
spatial trends in observed radial velocity and heliocentric distance, although we note that our 
goal in this work is not to model the Sagittarius stellar stream in precise detail. The massive 
dark matter component is disrupted on scales much larger than that of the systematic errors 
involved in both measuring and modeling the luminous debris, such that the general behavior 
of Sgr dark matter is grossly similar in both experiments, as we now address. 

The presence of a coherent stream of dark matter in the solar neighborhood is a robust 



prediction of both Sagittarius models, and would be a feature of any putative progenitor 
with cosmologically-motivated and observationally-consistent properties as discussed above. 
The material is stripped from the infalling satellite, and in both models the dark compo- 
nent of the leading tidal arm contributes a non-negligible fraction of dark matter to the 
solar neighborhood. We emphasize that the contribution of dark matter from Sgr is non- 
negligible despite the fact that the stellar arm of the tidal debris is displaced from the solar 
neighborhood by several kiloparsecs in the Milky Way disk plane in our baseline models. 
This occurs for two reasons. First, in any model in which the Sagittarius progenitor has a 
cosmologically-motivated progenitor halo, the progenitor halo is significantly larger than the 
stellar component of the progenitor galaxy and spawns a tidal stream that is significantly 
more spatially extended than the stellar debris stream. Indeed, the dark matter streams in 
our models are more than > 10 kpc in cross-sectional diameter. Second, the bulk of the 
dwarf's dark matter is stripped from the Sagittarius progenitor prior to the onset of stellar 
mass loss. The dark matter and stellar debris material do not follow precisely identical or- 
bits, so a cross-section of the dark matter stream at the plane of the disk is not necessarily 
concentric with a cross-section of the stellar debris. These features are both evident in Fig. 1, 
with the offset between the dark matter and stellar arms most prominent in the light Sgr 
model, since the progenitor in that case has been more heavily disrupted. 

Our two Sagittarius models shown in Fig. 1 could in principle underestimate the amount 
of Sgr dark matter near the Sun, because the axis of the leading stellar stream misses the 
solar neighborhood entirely in each case, with only the outskirts of the dark matter in that 
tidal arm contributing to a potential WIMP detection signal. In fact, previous models for the 
Sagittarius leading arm have indicated that it may impact the disk plane within a few kpc 
of the Sun [56], although observational explorations of the solar neighborhood have failed to 
detect any sub-population kinematically similar to streaming substructure with Sagittarius- 
like geometry within one or two kiloparsecs of the Sun [43, 44], and constraints on the 
descending portion of the leading arm are presently very poor (e.g., [55]; see also the "field 
of streams" depicted in ref. [41]). We note that the more recent investigation in ref. [57] of 
the Sgr stream in a triaxial halo suffers from a number of flaws, including the lack of a dark 
matter component in the progenitor as well as the adoption of a static Galactic potential. In 
contrast, our modeling techniques and initial conditions self-consistently resolve these issues 
in a numerical context that correctly treats the evolution of tidal disruption characteristics 
in the stellar debris. As a further point of distinction, we also emphasize that the preferred 
model of ref. [57] (in which the leading arm passes ~ 10 kpc from the Sun) requires a 
nearly-oblate Milky Way halo (at odds with the prolate shape preferred by long-lived warps 
in the atomic-hydrogen gas layer in the Milky Way [58] as well as the Holmberg effect of 
satellites being clustered around a plane normal to that of the central galaxy [59]) with 
an unusual orientation compared to the cosmological findings of refs. [60-62]. In any case, 
the orbital shape of the satellite's infall is sufficiently well-constrained that any debris model 
with a cosmologically-realistic dark subhalo will necessarily result in a wide stream of WIMPs 
raining into the solar system. 



2.2 Dark Matter Direct Detection Event Rates 

After obtaining the speed distributions characterized by f{v) in the left panel of Figure 2, 
for each of the four models we calculate differential event rates as a function of recoil energy 
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Figure 2. Left — Dark matter velocity distribution relative to Earth in the early Northern winter 
(day-number 7 relative to noon on December 31st, for the J2000.0 epoch), for each of our four models. 
The dotted histograms are best fit to a Maxwellian form by the thin solid lines in the upper panel, 
while the fractional differences from the Maxwellian fits are shown as a function of velocity in the 
lower panel. The excesses near v ^ 450 — 500 km/s in both of the Sagittarius models are real, 
and drift during the year from this wintertime peak to a summertime low of w ^ 400 km/s. Right 
— The integral quantity ,9(wmin) as a function of the minimum particle velocity, Wmin, required to 
initiate nuclear recoil in a direct detection experiment for each of our four models (upper panel; also 
shown is a comparison to the standard halo model (SHM) Maxwellian for a range of circular velocities 
Vc = 200 — 220 km/s) . The second panel from the top shows the residual of each model compared to 
a best-fit Maxwellian {not the SHM), the third panel from the top shows the residual of the models 
that include the Galactic disk compared to the host halo model, and the bottom panel shows the 
residual of the Sagittarius models compared to the halo + disk model. The integral g{vmin) is 
directly proportional to the event rate so the generally higher velocities in the models including the 
Galactic disk result in a factor of roughly two or more increase in event rate relative to the host halo 
model. The features induced by each Sgr stream model lead to further modifications of this prediction 
on the f^ 10 — 20% level {lower panel). For illustration, the alternate axes in the upper panel show 
corresponding nuclear recoil energies in unquenched keV,- for germanium and sodium target nuclei, 
assuming a WIMP mass of m^ = 10 GeV/c^. 



as in [2]: 
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where a^ is the WIMP cross-section for scattering on a proton (assuming here that the 
WIMP couples nearly equally to protons and neutrons), po is the WIMP density in the solar 
neighborhood, /i is the reduced mass of the proton and a WIMP with mass m^, A is the 
atomic mass number of the detector nuclei, F{E^) is the form factor of nuclear scattering as 
a function of recoil energy E^, and the quantity ^(fmin) is the integral in velocity space of 
the velocity distribution divided by the WIMP speed. 



5(^^min) = / ^dv. (2.2) 



Here, v^in = ^ / ^ i"' , where Ma and ha are the atomic mass and WIMP-nucleon reduced 

V >^A 

mass, and this is the minimum relative speed necessary for nuclear recoil to yield an energy E^ 
(lower WIMP masses require higher values of Umin at fixed recoil energy). In our analysis, we 
calculate g{vmin) using Earth rest-frame velocity distributions and the local WIMP density 
po directly from our simulations. We adopt the Helm form factor F{Et-) (as in [63]; see 
[64] for details and fitted parameters). Throughout, we choose an arbitrary value of a^ = 
10"^'' cm^ (= 10""^ pb) for absolute event rates. Predicted event rates can be scaled linearly 
for different values of a^-. We frame our results in terms of unquenched nuclear recoil energy 
E^, using the standard unit notation keVr, reminding the reader that the quenching factor 
for a particular detector material must be used to convert this to electron-equivalent recoil 
energy in keVee, i.e. -EkeVee = 9-^keVr- ^^^ ^^^ detector examples we investigate in this work, 
(g, x) ~ (0.199,1.12) for germanium [65], and {q,x) ~ (0.3,1.0) for sodium [11]. Note that 
we do not model the finite energy resolution of detectors. 

3 Results 

3.1 General Deviations from the Standard Halo Model 

Our primary results concern the non-trivial amount of dark matter donated to the solar 
neighborhood by the leading tidal arm of the Sagittarius dwarf galaxy (shown in Figure 1), 
but prior to this discussion, we briefly itemize pertinent features of the dark matter compo- 
nent of the primary halo emphasizing deviations from a Maxwellian distribution. Many of 
the gross features that characterize deviations from a Maxwellian have been pointed out in 
previous work [16-20]. We note that the particular systems we study have been tuned to 
equilibrium host halo or halo + disk configurations that represent many of the gross features 
of the Milky Way galaxy, but that these are not unique solutions for equilibrium models of 
the Milky Way. Dark matter direct search rates are usually interpreted in the context of a 
standard halo model (SHM) [14, 15] with p{r) oc r~^, a local density of po = 0.3 GeV/cm^, 
and a velocity distribution described by a Maxwellian form with a three-dimensional dis- 
persion (73D = \/3/2vc, where Vc is the circular speed at the solar radius which is assumed 
to be identical to the peak speed of the Maxwellian distribution. Fig. 2 shows that the 
equilibrium models we consider result in direct search rates significantly different from the 
specific Maxwellian assumed to derive SHM predictions (the shaded bands in the upper right 
panel). Furthermore, the velocity distributions deviate significantly not only from the SHM, 
but from the general form of a Maxwellian speed distribution. Fig. 2 shows comparisons be- 
tween velocity distributions in our equilibrium galaxy models and Maxwellian distributions 
with mean and dispersion that best fit the simulation data. The isolated host halo exhibits 
significant deviations from its best-fit Maxwellian distribution as shown in Fig. 2 and noted 
in previous studies and in agreement with the cosmological A'^-body results presented by 
ref. [18]. The same is true for the halo + disk and Sgr-infall models. In both cases best-fit 
Maxwellians underestimate the value of g(fmin) by ~ 20 — 40% over a wide range of recoil 
energies relevant to direct detection experiments. 

Relative to an isolated halo, the addition of the Galactic disk causes the host halo to 
contract globally as well vertically toward the disk mid-plane. The global contraction of the 
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Figure 3. Examples of predicted direct search event rates. We show time-averaged differential 
scattering event rate dR/dE as a function of recoil energy (keVr, unquenched), for scattering off of 
germanium {left panels) and sodium [right panels). In each detector, we show example spectra for 
light WIMPs of mass m^ = 5 and 10 GcV/c^. We show the fractional deviation in the event rate 
from the host halo model in the center panels, and the fractional deviation of the Sagittarius model 
event rates from the halo + disk values are shown in the lower panels. 



host halo in response to a disk shifts the velocity distribution to larger velocities (Fig. 2), 
resulting in relatively higher event rates (Figure 3). In addition to this gross shift, equilibrium 
Galaxy models including a disk exhibit speed distributions with significantly more broad, 
flat peaks. The kurtosis of the host halo speed distribution is iT ~ —0.3, while the galaxy 
models have K ~ —0.5 (by definition K = Q for a Gaussian distribution, and X ~ 0.1 for a 
Maxwellian). Scattering event rates are directly proportional to the local WIMP density po, 
which in all cases is higher for simulations including the Milky Way stellar disk than for the 
dark-matter-only host halo model. The solar neighborhood in this original halo resides in a 
region with WIMP density pQ ~ 0.53 GeV/cm^, while the three models including the stellar 
disk of the Milky Way have po ~ 0.61 — 0.63 GeV/cm^, a ~ 20% increase due to the overall 
contraction of the halo in response to the Galactic potential and the compression of the halo 
in the vertical direction due to the planarity of the disk. We note that our models have been 
tuned to Milky Way structural properties [38] (and also respond to variously-sized satellite 
impacts in a fashion consistent with that of the quiescent Galactic mass accretion history 
and the global spiral structure of the disk; see refs. [45, 66]), so these values of po illustrate 
the potential error latent in typical choices of this normalization. We present event rates 
calculated self-consistently with respect to the WIMP density po in each model; however, the 
changes in the speed distribution and, therefore, the integrated quantity (7(vmin)i represent 
far more significant alterations to predicted event rates. Figure 2 and Figure 3 show that the 
shift in the speed distribution increases ^(fmin), and thus event rates, by a factor of several at 
high Umin; a boost that is particularly important for light WIMP masses m^ ~ 5 — 20 GeV/c^ 
suggested by several recent experiments [10-12]. 

Unsurprisingly, the gross effect of the Galactic component is to deepen the gravitational 
potential, resulting in generally larger relative speeds of dark matter particles. It is interesting 



to assess whether this offset can be modeled simply so that they may be incorporated into 
future analyses, such as those proposed by ref. [67] among others. A simple proposal would 
be to employ a model for adiabatic halo contraction [68, 69] on the original host halo, which 
is well described by the standard NFW profile form. Specifying a velocity distribution would 
still be a challenge, but a simple proposal would be to employ the standard Eddington relation 
on the contracted halo (see the relevant exercises in ref. [70] , for example) . Such an approach 
would certainly not be self-consistent as both formalisms assume spherical symmetry and 
the Eddington relation does not yield a unique speed distribution; however, it is interesting 
to explore such an option as a simple, practical alternative to perform a gross correction 
to account for contraction without the expense of constructing equilibrium models of the 
Galaxy and its halo.^ We find that our models can be described in this way at a similar level 
of precision as using the best-fit Maxwellian for each model. In particular, we find residuals 
similar to the halo + disk residuals for the best-fit Maxwellian in Figure 2, with adiabatic 
contraction modeling resulting in residuals from ~ 20 — 50% for fmin > 300 km/s. At this 
level of precision, a simple analytical correction may be able to map results from a host halo 
identified in a dark-matter-only simulation to event rates in the solar neighborhood. An 
even more parsimonious exercise recovers the halo/galaxy speed distribution by scaling each 
particle's velocity in the host halo by an amount equal to the increase in the distribution's 
peak speed; for our models, this simple adjustment produces event rates within ~ 5% of 
those yielded by the halo + disk analysis. 

3.2 Dark Matter from Sagittarius at Earth 

Our most novel results pertain to the effects of the Sagittarius debris on direct search rates. 
Before exploring the influence of the Sagittarius tidal debris in detail, several points are 
worthy of reiterating here. Our Sagittarius models are cosmologically- motivated and are 
designed to bracket the range of halo masses that a Sgr-like galaxy would be expected to have 
as it merged with the Milky Way. Both Sagittarius models produce debris streams that are in 
broad agreement with the known morphological and kinematical properties of the observed 
Sagittarius stream stars. Moreover, as noted already in Fig. 1, in these self-consistent models 
of Sagittarius evolution, the dark matter stream is signiflcantly wider spatially than the 
stellar stream, spanning > 10 kpc in the direction transverse to the stream. Lastly, the dark 
matter and stellar debris streams are not necessarily coincident; in fact, at the point where 
the Sagittarius stream penetrates the disk closest to the solar neighborhood, the dark matter 
and stellar streams are not coaxial. The offset between the dark matter and stellar streams 
stems from the fact that the progenitor Sagittarius dark matter halo is significantly more 
extended than the stellar component of the progenitor dwarf galaxy, so that dark matter 
is typically liberated prior to stellar material. Consequently, the dark matter and stellar 
streams do not follow the same orbits. These considerations suggest that the dark matter 
component of the Sagittarius stream may likely be a very important contributor to scattering 
in earthbound direct detection experiments even if the stellar stream penetrates the Milky 
Way disk more than two kpc away from the Sun as is currently suspected (for example, in 
the modeling of ref. [57], which incorporates the general observational results of [42, 55], 
among others). These broad points represent an important addition to the literature on the 
influence of Sagittarius for dark matter direct detection experiments. 



^Among other caveats to such a formahsm, some numerical experiments of Galaxy-scale halos have indi- 
cated that halos may expand as baryonic feedback ejects material from the galaxy [71, 72], working against 
the tendency toward adiabatic contraction in a way not easily modeled by analytic forms. 
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Figure 4. Annual modulation of the dark matter scattering event rate. Upper left panel: at a fixed 
value of Wmin — 430 km/s, corresponding to a recoil energy of E,- ~ 4.6 keVr for a germanium detector, 
assuming a light- WIMP mass of m^ — 10 GeV/c^, we show the best-fit amplitude A and peak day- 
number for each model. Upper right panel: The fractional annual modulation amplitude as a function 
ofwmin, i-G. (max — min)/(max -|- min) . Lower panels: The dependence of the peak day- number on the 
minimum velocity Umin, with equivalent recoil energies for germanium- and sodium-based detectors 
plotted on alternate upper axes for m^ — 10 GcV/c^. Note that the peak day- number shifts phase 
to Northern winter, for small unquenched recoil energies E,- < 1 (2) keVr for germanium (sodium), 
corresponding to Wmin — 0.89wc in agreement with ref. [73]. In the lower right panel, the shaded regions 
represent the ±lcr ranges in the modulation peaks determined by DAMA/LIBRA for scattering on 
sodium, i.e. peak day-numbers of 136, 142, 146 ± 7 for recoil energy bins E^ = 2 — 4, 5, 6 keVee [11]. 
Where present, thin colored lines represent the maximum possible signal induced by each Sgr model, 
as discussed in the text. 



The overall fraction of dark matter in the solar neighborhood is relatively small in all 
of our models, varying from ~ 1 — 2% of po and may reasonably be as high as ~ 5% of po 
if po = 0.3 GeV/cm^ rather than the somewhat higher values in our models. Consequently, 
integrated event rates are altered by only relatively small amounts by the stream. Figure 2 
and Figure 3 show that the stream debris can alter event rates by up to ~ 10 — 20% at 
relative speeds greater than the typical relative stream speed, fmin ^ 400 km/s (if the local 
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density from the host halo is as low as /jq = 0.3 GeV/cm^ as in the SHM, this boost can 
be as large as ~ 25%). Note that the Sagittarius debris give a relative enhancement to the 
speed distribution at speeds near v ~ 400 — 500 km/s in Fig. 2. The best-fit Maxwellian is a 
better description of the models that include Sagittarius debris than it is of the halo + disk 
models because the high-speed particles in light and heavy Sgr influence the fitted kurtosis 
sufficiently that some of the halo's natural non-Maxwellianity is accounted for more closely. 

Rather than representing large modifications to the overall annually-averaged event rate, 
the speed distribution features around 400 — 500 km/s due to the Sagittarius debris represent 
potentially interesting peculiarities that may be explored with the annual modulation of the 
event rate by current-generation experiments as well as directional detection efforts address- 
ing the North Galactic Pole (latitude h ~ +90° with respect to the solar neighborhood) in 
the long-term future. The stripped dark matter from Sgr falls coherently toward the solar 
system from the North Galactic Pole at a speed of order ~ 400 — 500 km/s. The speed 
relative to the Earth peaks during Northern winter, nearly opposite in phase to the relative 
speed between the Earth and the dark matter in the primary host halo of the Milky Way, 
due to the geometry of the Sgr orbit with respect to the Milky Way disk plane. 

We present results regarding the annual modulation of dark matter direct search rates in 
Figure 4. In the upper left panel of Fig. 4, we show an example of the direct search event rate 
as a function of day-number (relative to noon on December 3f st during the J2000.0 epoch) for 
a germanium detector, a light WIMP with mass m^ = 10 GeV/c^, and unquenched nuclear 
recoil energy of ~ 4.6 keVr- The amplitude of the annual modulation of (?(fmin) is shown 
in the upper right panel of Fig. 4. We show the amplitude as the difference between the 
maximum and minimum event rates achieved within a yearly cycle divided by the sum of the 
maximum and minimum rates. If the modulation induces a simple sinusoid superimposed 
upon a constant background, this yields the amplitude of the sinusoid. The amplitude of the 
annual modulation signal is generically an increasing function of f min because the number of 
particles in the tail of the speed distribution is very sensitive to small shifts in the central 
position of the speed distribution. The SHM predicts a much more shallow rise in the 
fractional modulation amplitude at relatively high speeds than that yielded by all four of our 
models, as shown in the upper right panel of Fig. 4. Each of the models containing a stellar 
disk results in a fractional amplitude that is more than 20% lower than that of the host halo 
model alone. 

The Sagittarius debris alters predictions for the amplitude of the annual modulation 
markedly. In particular, both of our fiducial models of Sagittarius {dashed lines) exhibit an 
annual modulation amplitude as much as ~ 20 — 30% less than the modulation amplitude in 
the Galaxy models with no Sagittarius debris (and a factor of ~ 2 smaller than in the host 
halo) in the range 400 km/s < Wmin ^ 500 km/s. The modulation amplitude is reduced in 
the Sagittarius models because the dark matter stream of Sagittarius rains down upon the 
Galactic disk from the North Galactic direction, so that this component is distinctly out of 
phase with the background "wind" of dark matter particles from the primary Galactic halo 
[21-24]. 

As we have alluded to previously, the location of the impact of the leading stellar stream 
of the Sagittarius dwarf on the Milky Way disk relative to the solar neighborhood remains 
uncertain, partly because high-latitude data must be extrapolated to the disk plane in order 
to estimate the impact position [55], and partly due to uncertainty in the Sun's position. 
Consequently, we have estimated the possible range of annual modulation effects that may 
be attributable to the Sagittarius debris by artificially shifting the position of the Sagittarius 
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impact with respect to the Sun in our models by amounts consistent with observational 
constraints on the relative position of the Sagittarius impact with the disk and the Sun. 
As shown in the lower panels of Figure 1, the solar neighborhood in the light Sgr model 
is serendipitously near the peak of the Sagittarius dark matter density (despite the stellar 
stream being several kpc away, because the central axes of the dark and luminous tidal arms 
are not coincident), while the heavy Sgr model's solar position is ~ 2 — 3 kpc farther from 
the Galactic Center than the axis of the dark stream. Our models suggest that if the Sun 
is ^ 10 kpc or more from the center of the stellar stream impact on the Galactic plane, the 
influence of Sagittarius on direct search experiments may be quite small; however, current 
modeling and observational constraints indicate that the leading arm is probably significantly 
closer to the solar neighborhood [55]. 

To probe this uncertainty, we estimate a range of possible Sgr-related effects by taking 
the maximal local dark matter density commensurate with observational determinations of 
the stellar stream position. To a good approximation, the maximal influence of Sagittarius 
occurs when the leading dark matter arm falls directly onto the solar position. This can 
be achieved in both models even when the peak density of the stellar arm is more than a 
kiloparsec from the Sun (see Fig. 1). Moreover, we note that although observations indicate 
that the solar neighborhood is not significantly contaminated by former Sagittarius stars 
[43, 44], kinematically-streaming sub-populations that may be associated with Sagittarius 
would only contribute on the order of ~ 1% to the local stellar density, making them very 
difficult to rule out on heliocentric-distance scales larger than a kiloparsec. We refer to the 
cases in which the relative position of the Sun maximizes the annual modulation influence of 
Sagittarius as our maximal models for the heavy Sgr and light Sgr simulations respectively. 

The maximal Sagittarius cases are shown by thin, solid lines in the upper right panel 
of Fig. 4. The additional decrease in the annual modulation amplitude is significant. In 
particular, the annual modulation amplitude can be decreased by as much as a factor of two 
in the maximal Sagittarius models compared to the halo + disk models in the speed range 
400 km/s < ^min < 500 km/s. 

The Sagittarius stream dark matter has an important influence on the phase of the 
annual modulation of scattering rates. The phase of the annual modulation as a function of 
f^mirn specified by the day of the peak signal (with day-number 1 set to noon on December 31st 
in the J2000.0 epoch), is shown in both lower panels in Fig. 4. The peak days in the host halo 
and halo + disk models are close to the canonical value of day 152.5 for t^min ^ 190 km/s. 
At lower values of Vmiru the phase switches to a peak day in the Northern winter, close 
to day 335 (or equivalently day —30 as shown in Fig. 4). The shift in the phase of the 
oscillation at low energies occurs because the shift to higher relative speeds caused by the 
greater relative speed of the Earth compared to the primary host halo dark matter particles 
during Northern summer augments the high-speed portion of the speed distribution, but 
depletes the low-speed portion of the speed distribution so that low-energy scattering rates 
are reduced during Northern summer. Conversely, low-energy scattering rates are enhanced 
during Northern winter. For the SHM, the peak day at high v^\^ occurs on day 152.5 (June 
2nd), while the peak day at low Vmin occurs on day 335 (December 1st). The shift of the 
peak day from Northern winter to Northern summer occurs at v^am — O.SOvc in the SHM 
[73]. All of our models show a similar phase reversal in the peak day, though the value of 
Vmin at which the shift occurs varies slightly with respect to the SHM value, as we discuss 
shortly in the context of the Sagittarius models. 

The Sagittarius streams in our simulations alter the annual modulation signals in non- 
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Figure 5. The phase reversal of the annual modulation signal at relatively low recoil energies, 
for each of our four models. Note that the peak occurs in Northern winter for fmin < 0.89wc as 
expected[73], for the host halo and halo + disk models. The Sagittarius models reverse phase at lower 
I'min- This discrepancy could potentially manifest itself following a large number of scattering events 
observed by the next generation of direct detection experiments, which may feature sub-keV recoil 
energy thresholds and improved energy resolution. As in other figures, equivalent recoil energies for 
germanium- and sodium-based detectors are plotted on alternate upper axes for m^ = 10 GeV/c^. 
Heavier WIMPs, near the canonical mass range of m^ ^ 10^ — 10^ GeV, imply higher recoil energies 
than those shown on the upper axes here. 



negligible ways. The differences between the halo + disk or host halo models as compared 
to the Sagittarius models depend upon fmm- The most significant deviations are for Vmm ^ 
250 km/s and the range of speeds typical of the relative speed of the Earth with respect 
to the Sagittarius stream, 400 km/s < Vmm ^ 500 km/s. In the higher speed range, the 
maximal shifts in peak day bring the peak between 5 and 12 days earlier in the year for the 
fiducial heavy Sgr and light Sgr models respectively. In the maximal models, the shift can 
be larger, bringing the peak day as much as 12 and 18 days earlier in narrow energy ranges 
in the heavy Sgr and light Sgr models respectively. The shifts in peak day are comparable 
in magnitude to the errors on the DAMA/LIBRA measurement of annual modulation [11]. 
Moreover, the shifts in our models bring the peak day earlier by amounts and with energy 
dependence comparable to the offset between the SHM peak day and the peak day quoted by 
the DAMA/LIBRA experiment, for which the signal peaks at day-number 136, 142, 146 it 7 in 
the recoil-energy bins 2 — 4, 5, 6 keVee, as indicated by the shaded regions in the lower right 
panel of Figure 4. Additionally, the energy dependence of the annual modulation amplitude 
in our models may partially mitigate the discrepancy between the low amplitude observed by 
DAMA/LIBRA at energies 2-6 keVee {A ~ 2% of the mean rate 0.0116 cpd/kg/keVee) and 
the somewhat higher amplitude A ~ 13% quoted by CoGeNT at energies 0.5 — 3.0 keVee. 
Conversely, the peak day shift induced by the Sgr debris becomes indistinguishable from 
noise, if the solar position can be located as much as ~ 5 kpc from the center of the stellar 
debris arm, a distance somewhat greater than contemporary estimates. 

It is tempting to suggest that experiments such as DAMA/LIBRA or CoGeNT may be 
probing the stream of the Sagittarius galaxy, but any such statement is subject to numerous 
important caveats. First, the DAMA/LIBRA results are controversial and are challenged by 
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other experiments. Second, the DAMA/LIBRA peak day estimates deviate from the SHM 
(and halo + disk model) predictions by only 2a or less. However, another caveat to any such 
statement strikes at a limitation of our simulation analysis. As we have already mentioned, 
the position of the Sagittarius stream impact on the disk remains uncertain. Further the 
angle at which the Sagittarius stream impacts the disk is also uncertain. Finally, despite 
the fact that these simulations are the most cosmologically complete of such efforts, having 
accounted for realistic dark matter halos in each galaxy, some small discrepancies remain 
between our results and observations of the sky-position and radial velocity of the stellar 
debris. Unfortunately, a parameter search using numerical simulations is not yet practicable. 
As an example of the concrete effects of this uncertainty on our predictions, the amount of the 
phase shift induced by the Sagittarius stream and, in fact, whether the shift is toward earlier 
in Northern spring (as shown in Fig. 4) or, perhaps, later in Northern summer is sensitive 
to the angle at which Sagittarius impacts the disk. Contemporary data and techniques do 
not suffice to specify this impact direction, so it is not possible to interpret our results as 
a firm statement that the Sagittarius stream causes the peak day to occur earlier by some 
precise amount (the amplitude of the annual modulation is not subject to this caveat). A 
proper interpretation of our simulation results is that we have shown that it is possible, in 
self-consistent models designed to mimic the Milky Way with the Sagittarius impact, for the 
Sagittarius debris to induce significant, energy-dependent shifts in the annual modulation 
phase even when the stellar stream of the Sagittarius debris remains several kiloparsecs from 
the Sun. 

We now turn to the phase reversal at low velocities/recoil energies. As we have already 
mentioned, the phase reversal of the annual modulation at low energies is a well-understood 
effect, occurring at Vmm — 0.89vc in the SHM [73] and at similar values in our host halo and 
halo + disk models. As depicted in Figure 5, the presence of Sgr dark matter in the solar 
neighborhood affects not only the phase of the overall modulation, but also the value of Umin) 
and thus the recoil energy, at which the phase reversal occurs. In particular, the Sagittarius 
models reverse phase at speeds ~ 10 — 15 km/s lower than the host halo and halo + disk 
models. For relatively low-mass WIMPs with rriy. = 10 GeV, this corresponds to a shift in 
recoil energy at the phase reversal point of approximately Aii^r ~ 0.1 — 0.2 keVr. For larger 
WIMP masses in the regime m^ > 100 GeV/, this shift in the reversal energy would be larger 
by a factor of 10 or more. 

This Sagittarius-induced feature may be probed by future detectors with low energy 
threshold and improved energy resolutions and may be one of the distinguishing features 
of Sagittarius if WIMP astronomy can ever be undertaken in an era with very large direct 
search event rates (as explored by ref. [67]). Exploring this particular signature is well 
suited for efforts to develop very low-threshold direct search detectors with greatly improved 
energy resolution, examples of which include extensions of germanium-based experiments like 
MAJORANA [74, 75], advanced threshold-lowering and background-eliminating technologies 
like those proposed by the CDEX-TEXONO and SuperCDMS/CDMSLite cohaborations (as 
in refs. [76-78], respectively). This avenue of exploration may be particularly fruitful in 
the future as the local circular speed Vc of the Milky Way is refined by next-generation 
astrometric surveys so that a similarly tight constraint can be placed on the phase-reversal 
energy predicted by the SHM and similar models without Sagittarius debris. 

There is a final addendum to our model results that is worth stating explicitly. In our 
models, the potential interior to the solar position is dominated by stellar mass, rather than 
dark matter. In the event that the local dark matter density contributed by the primary host 
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halo of the Milky Way (as opposed to the Sagittarius stream) is as low as the SHM value 
of po ~ 0.3 GeV/cm^, we expect viable models of Sagittarius evolution to remain broadly 
similar. The implication is that the relative influence of Sagittarius on direct search scattering 
rates could be yet larger than we have estimated here. Although we have not performed a self- 
consistent model that results in such parameters for the Milky Way halo, it is interesting to 
comment on how our results would change in such a scenario. If the local density contributed 
by the host halo of the Milky Way were as low as po ~ 0.3 GeV/cm^ (rather than roughly 
twice this value as in our equilibrium halo/galaxy models), the Sagittarius debris could 
contribute as much as ~ 5% of the local dark matter density and the direct search signatures 
we have explored would change as follows: the phase-shift in the annual modulation signal 
could be as large as ~ 20 — 25 days, compared to the ~ 10-day shift yielded by our fiducial 
models, and the fractional amplitude of the annual modulation signal could remain as low as 
~ 5% even at relatively high speeds near Umin ~ 500 km/s (well below the SHM prediction 
as well as each of the cases we examine here). In addition, the change in the recoil energy at 
which the modulation undergoes phase reversal could be as much as twice the fiducial change 
we show in Fig. 5. These estimates, based on a scenario in which the solar neighborhood 
is near the peak of the dark matter stream and the local density contributed by the parent 
halo is as low as po = 0.3 GeV/cm^, likely represent maximum plausible influences that 
Sagittarius stream material could have on direct search rates without fine-tuning. 

4 Discussion 

We have studied predictions for dark matter direct search scattering rates within the context 
of isolated numerical models of a Milky Way-like system designed to reproduce the basic 
properties of the Galaxy, including models of the infall, merger, and tidal disruption of 
the Sagittarius dwarf system. In modeling such specific features, isolated simulations of 
this kind complement large-scale cosmological simulations. In accord with previous high- 
resolution studies of cosmological dark matter halos [16, 18-20], we find that deviations from 
standard halo model (SHM) assumptions in observationally viable model Milky Way systems 
can significantly alter direct search rates relative to SHM predictions. In agreement with 
these studies, equilibrated host halo systems and halo systems in equilibrium with a Galaxy 
exhibit speed distributions that differ markedly from the Maxwellian form, in particular being 
significantly platykurtic. 

Not surprisingly, in all three models that include an equilibrated Galactic stellar disk, 
we find that the increased relative speeds of dark matter particles caused by the additional 
acceleration provided by the disk and the contracted halo result in significantly larger scat- 
tering rates. The precise value of this enhancement can be large (a factor of several) and is 
energy dependent (Fig. 2 and Fig. 3). Studies based upon A^-body realizations of Galaxy- 
analog halos in cosmological simulations make the implicit assumption that a Galaxy-sized 
halo in a cosmological numerical experiment involving only dark matter will faithfully re- 
flect the solar neighborhood in the real Galaxy [19, 20]. For practical purposes, mapping 
cosmological N-body results onto equilibrium models containing a galaxy is non-trivial. We 
have shown that our equilibrium halo + disk models can be used to predict rates to within 
50% by contracting the dark matter halo using standard adiabatic contraction techniques 
[68-70]. Meanwhile, simply scaling the speed distribution of the A^-body host halo by the 
mean velocity can reproduce the halo/galaxy distribution to within ~ 5% precision, signi- 
fying that dark matter-only predictions for speed distributions can be mapped to models 
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that include a Milky Way Galaxy by scaling speeds up to the rotation speed in the solar 
neighborhood. These results are broadly commensurable with previous A^-body work and a 
significant caveat to these results is that our equilibrium models represent only one possible 
equilibrium solution for the halo/galaxy system that is not unique and does not result from 
self-consistent cosmological evolution. 

Our most novel results pertain to the influence of Sagittarius debris material on pre- 
dicted direct search event rates. Our models demonstrate that the Sagittarius stream debris 
can have an important influence on direct search scattering rates even when the stellar stream 
of the Sagittarius debris is centered several kiloparsecs from the solar neighborhood, as is 
thought most likely based on contemporary analyses [42, 55, 57]. The reasons for this are 
twofold. First, the Sagittarius dark matter stream is significantly broader than the stellar 
stream in our models. The Sagittarius stream is a non-negligible contribution to the nearby 
dark matter content over an area many kiloparsecs in diameter in the plane of the Galactic 
disk (Fig. 1). Second, the Sagittarius stellar and dark matter streams are not spatially co- 
incident in our models, having drifted away from co-axiality during evolution in the Milky 
Way's tidal field. The peak of the dark matter density contributed by the stream impacts the 
Galactic disk several kiloparsecs from the peak of the stellar density, and the spread in sur- 
face density demonstrated by Fig. 1 indicates that expected event-rate boosts for detection 
experiments should be important out to this distance, in our adopted formalism. 

We reiterate an important caveat before continuing; although the spatial and kinematic 
distributions of Sagittarius stellar debris in our models are generally good matches to the 
observational properties of that debris, as elaborated in [45], some discrepancies do remain. 
Moreover, the relative position of the solar neighborhood with respect to the Sagittarius 
stream's impact point on the Galactic disk is still poorly constrained. The simulations we 
analyze are among the most complete descriptions of the Sagittarius debris in the literature 
and an exhaustive search of the initial parameter space is not computationally feasible, so 
these uncertainties cannot be explored in detail. A proper interpretation of our results would 
be that the Sagittarius dark matter debris may give rise to significant signals in direct-search 
experiments even if the Sagittarius stellar debris is confined to a distance of several kiloparsecs 
from the solar neighborhood. Tuning the stellar stream more finely would not impact our 
results, insofar as variation in Sgr WIMP surface-density at the Sun is less important to 
the event-rate calculation than is the vertical velocity distribution of those WIMPs. The 
important and general implication of this result is that the stellar stream of Sagittarius 
being several kiloparsecs from the solar neighborhood does not preclude Sagittarius from 
having a significant effect on dark matter experiments. Furthermore, near-future constraints 
on the location of the stellar debris may not suffice to preclude significant dark matter from 
Sagittarius in the solar neighborhood. 

The Sagittarius stream gives rise to several important effects on dark matter experi- 
ments. The high relative velocity of the Sagittarius dark matter stream relative to the Earth 
boosts the rate of high-energy recoil events by ~ 20% in our models and perhaps as much as 
~ 40 — 45% depending upon the local dark matter density contributed by the primary halo of 
the Milky Way (Fig. 2 &; Fig. 3). Sagittarius reduces the annual modulation amplitude by an 
energy-dependent factor that may be as large as a factor of two at fmin ~ 420 km/s (Fig. 4) 
relative to models with no Sagittarius debris. This energy-dependent suppression could help 
explain the disparate values found by DAMA/LIBRA (where the modulation amplitude is 
A ~ 2% of the mean low-energy rate dR/dE ~ 0.0116 cpd/kg/keVee in the energy bin 
2 — 6 keVee = 6.7 — 20 keVr) and CoGeNT (having an amplitude A ~ 13% of the mean rate 
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for recoil energies ~ 0.5 — 3.0 keVee = 2.3 — 11.3 keVr) as well as their discrepancies compared 
to the amplitudes expected within the SHM [11, 12], because our Sagittarius models result in 
amplitudes that increase much more sharply with recoil energy in this range than the SHM 
formalism as well as the host halo and halo + disk models. 

The geometry of the Sagittarius impact on the solar system causes the signal from the 
Sagittarius stream to peak during Northern winter, in agreement with previous studies [21- 
25] (as well as the general study of debris flows in ref. [79]). Our models of Sagittarius yield 
recoil energy-dependent shifts in peak day number of between ~ 5 and 25 days earlier in 
the year than the SHM peak day-number of 152.5. Both DAMA/LIBRA and CoGeNT have 
indicated a similar behavior, with both experiments finding trends between peak day-number 
and recoil energy [11, 12]; however, the observational error in the peak remains on the order of 
a few days and the uncertainty in simulation programs that aim to model Sagittarius remain 
significant. Nevertheless, our models suggest such a shift is reasonable given contemporary 
knowledge of Sagittarius debris structure. The Sagittarius streams in our models also cause 
the phase reversal energy to be lowered (Fig. 5). Exploiting this signature, in particular, 
to help identify dark matter or use dark matter searches to perform WIMP astronomy will 
benefit greatly from future low-threshold detectors with improved energy resolution, such as 
those being considered by refs. [76-78]. 

There have been a number of previous studies of the influence of Sagittarius debris on 
direct dark matter searches, including refs. [21-26]. By and large, these studies have utilized 
the SHM formalism and contrived approximations for the local distribution of Sagittarius 
debris based on observational limits on the Sagittarius stellar debris contribution to the solar 
neighborhood. Further, the prospects for Sagittarius to significantly influence direct search 
experiments have been challenged as contemporary data suggest that the stellar component 
of the Sagittarius debris is centered a few kpc from the solar neighborhood [55, 57] (and here 
we reiterate that although refs. [43, 44] find no evidence of large coherent sub-populations 
in the solar vicinity, percent-level debris flows are presently unconstrained beyond one or 
two kiloparsecs from the Sun). Our study complements previous work in many respects. 
Most importantly, we have analyzed a self-consistent simulation of Sagittarius dwarf galaxy 
accretion involving a realistic dark matter component, doing so in a controlled and isolated 
simulation such that event-rate implications can be clearly discerned. Crucially, we have 
shown that the proximity of the Sun to the stellar stream alone cannot necessarily be used 
as an indicator of the local dark matter contribution from Sagittarius, as the dark matter 
flow accompanying the luminous debris is much more widely spread across the Milky Way 
disk and not necessarily coincident with the stellar material. 

Comparing to previous work in more detail, we present the flrst analysis of the influ- 
ence of Sagittarius dark matter based on self-consistent models of the Sagittarius infall that 
describe the observed Sagittarius debris within observational uncertainty. Unlike previous 
studies, we emphasize that the Sagittarius debris induces an increase in the event rate by 
~ 10 — 20% in our fiducial models (as much as ~ 40% in our maximal models), and that 
annual modulation fractional amplitudes are diminished by ~ 20% — 50% in the presence of 
that debris at Earth. The phase shifts in the annual modulation that we find are somewhat 
larger than those presented by refs. [21, 22] (for comparable dark matter density contribu- 
tions), due to the relatively higher speed attained by Sgr particles in our modeling. We 
find qualitatively similar phase behaviors to those of refs. [23, 24], near the reversal point 
at which the modulation amplitude changes sign. Generally, our work agrees well with past 
commentary on the detectability of WIMP streams in current- and next-generation detec- 
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tion efforts [25-27], and specifically we identify the Sagittarius dark matter stream as an 
achievable target for direct-search science over the next decade. 

The effects of Sagittarius that we describe in this manuscript may be relevant to dark 
matter searches generally. However, if the dark matter is indeed relatively light {m^ < 
20 GeV as we have assumed in our illustrative examples), the effects of Sagittarius debris 
on scattering rates are particularly important because of the large relative speed of the 
debris stream at the Earth. In either case, future direct search experiments may probe such 
signatures, though a future generation of low-threshold detectors with fine energy resolution 
[76-78, 80] may be necessary in the event that the dark matter mass falls in this lower 
range. In either case, our analysis suggests that the effects of Sagittarius debris on direct 
search experiments will not be negligible given contemporary limits on the position of the 
Sagittarius stellar stream. In the far future, the features induced by Sagittarius debris may 
be among the early measurements to be made in an era of WIMP astronomy with large direct 
search rates [67]. 
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